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Abstract. Recent studies have shown that for suitable initial conditions both super- and sub- 
Chandrasekhar mass carbon-oxygen white dwarf mergers produce explosions similar to observed 
SNe la. The question remains, however, how much fine tuning is necessary to produce these 
conditions. We performed a large set of SPH merger simulations, sweeping the possible parameter 
space. We find trends for merger remnant properties, and discuss how our results affect the 
viability of our recently proposed sub-Chandrasekhar merger channel for SNe la. 
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1. Introduction 

The nature of type la supernova (SN la) progenitors, and their evolution leading up 
to the explosion, are still not yet well understood ( Van Kerkwijk et al.||2010 ). The most 
widely accepted scenario is that of a carbon-oxygen white dwarf (CO WD) in a binary 
that acquires mass through either accretion from a non-degenerate companion (the single- 
degenerate channel) or merger with another WD (the double-degenerate channel), until 
it approaches the Chandrasekhar mass, M ch ( |Howell|201lT ). This scenario still has several 
unresolved issues, including discrepancies between the rates of SNe la and progenitors 
( Maoz 2008 Ruiter et al. 2009 1 , the lack of sufficient numbers of supersoft sources for 
the single-degenerate (and potentially double-degenerate) channel ( Di Stefano||2010a|b| 



Gilfanov fc Bogdan||2010[ ) , and the still-poorly understood nature of the deflagration-to- 



detonation transition, required for Mch explosions to look like observed SNe la (Howell 



2011 1. These issues suggest the utility of considering alternate SN la channels, for instance 



the proposed channels of Pakmor et al. (20101, [Guillochon et al. (20101, and, for this 



paper, 



Van Kerkwijk et al. (2010). 



| Van Kerkwijk et al. propose an SN la channel that involves the merger of two CO 
WDs whose total mass is below Mch- As shown in Loren-Aguilar et al. ( 2009 ), the merger 
between two equal mass 0.6 M Q WDs results in a massive, uniformly rotating core with a 
differentially rotating, sub-Keplerian, fat disk. While the merger is hottest at the centre of 
the core, it is not hot enough to ignite runaway carbon fusion. The disk then accretes onto 
the core on a characteristic timescale of several hours. The accretion compresses the core, 
driving up temperatures until the accretion and carbon fusion timescales are comparable. 
The neutrino cooling timescale is six orders of magnitude longer, and therefore a carbon 
fusion runaway results. 
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This channel possesses a number of advantages. If most merging CO WD binaries 
(regardless of total mass or mass ratio) are assumed to result in SNe la, it would account 
for the observed SN la rate (Van Kerkwijk et al. 2010). Since lower-mass CO WDs 



generally come from lower-mass stars, it also explains the observed decrease of typical 



SN la luminosity with increasing host stellar population age. Sim et al. (20101 recently 



showed that sub-Chandrasekhar detonations can produce light curves and spectra very 
similar to observed SNe la, without resorting to a deflagration-to-detonation transition. 
This channel stands in contrast to He-CO sub-Mch merger channels, which involve the 
detonation of a helium layer around a sub-Mch CO WD, producing light curves and 
spectra that do not resemble observed SNe la (Sim et al. 2010 |Howell|[2011[ see e.g. 
Woosley fc Kasen||2011 on whether or not the effect of helium can be minimized) . 



In order to determine the viability of | Van Kerkwijk et al. s analytical sketch, it must be 
determined which regions of the CO WD merger parameter space will lead to detonations, 
and whether, e.g., nearly equal mass WDs are required, and, if so, how close to unity the 
mass ratio must be. 

We have run a series of 45 high-resolution smoothed-particle hydrodynamic (SPH) 
simulations of white dwarf mergers, exploring a broad range of parameter space to seek 
global trends in merger remnant properties as a means of interpolating between numerical 
results. Our results can be used to create initial conditions for models of the further 
viscous and thermal evolution of the remnants (e.g. Shen et al.|2011 ) to determine if any 
of them reach conditions necessary for carbon detonation. They also have utility outside 
our own project, as they would allow anyone to estimate the properties of a CO WD 
merger remnant without having to simulate the merger themselves. 

This report focuses on the implications of our results for the sub-Mch merger channel. 
Details regarding trends and a possible semi-analytical picture of CO WD mergers will 
be presented in a forthcoming publication. 



2. Code and Initial Conditions 

We used the SPH code Gasoline ( Wadsley et al.(2004 |. Gasoline uses the formulation 
for artificial viscosity in Monaghan ( 1992 ), augmented with a |Balsara ( 1995 ) switch and a 
method that reduces the artificial viscosity coefficient a when no shocks are present. We 
have modified Gasoline to use the Hclmholtz equation of state ( Timmes]|2010 details in 
Timmes fc Swesty||2000 1 . No nuclear reaction network was used, as previous simulations 



have shown that nuclear processing is not important to the hydrodynamic evolution of 
lower-mass mergers ( |Guerrero et al.|2004| [Toren-Aguilar et al.|2009[ ). A n in situ detona- 
tion, however, may arise during mergers between nearly equal-mass CO WDs where the 
accretor has mass > 0.9 M Q ( jPakmor et~alT1|2010[ [20TT] ) . 

Our input binaries are composed of 8 x 10 4 to 4 x 10 5 particles, depending on their 
mass. White dwarfs were relaxed individually and given equilibrium temperatures of 
~ 10 7 K. Their initial separation is chosen such that the lower-mass binary fills its Roche 
lobe. Since the synchronization timescale is expected to be much longer than the angular 
momentum loss timescale due to gravitational radiation, our binary is unsynchronized at 
the start of merger (Marsh et al.||2004 but see their discussion and Fuller fc Lai||2011). 



3. Results 

The merger remnants for all our runs consist of a core and a sub-Keplerian disk. For 
highly unequal mass mergers, almost the entire accretor forms 85 - 90 % of the core, while 
~ 70% of the donor forms the disk; for equal mass mergers, < 10% of donor and accretor 
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Figure 1. Merger remnant maximum temperature and corresponding density. Black points are 
merger remnants, and lines connect points with the same accretor mass (from bottom to top, 
0.4, 0.5, 0.55, 0.6, 0.65, 0.7, 0.8, 0.9 and 1.0 Mq). Maximum temperatures of nearly equal mass 
mergers have been adjusted to account for mixing in their convectively unstable cores. Also 
shown are contours of neutrino cooling timescale t v , carbon fusion timescale r cc and entropy S. 
The t v = t cc and Td yn = r cc lines denote where the carbon fusion timescale balances the neutrino 
cooling and dynamical timescales, respectively. The P = 2P(T=0) line marks the upper bound 
to where matter is degenerate. 



go into the disk. Other mergers fall on a continuum between these two extremes. Some 
results important to the possibility of the remnant exploding are summarized below: 

• The maximum temperature T max of the system is nearly independent of mass 
ratio, and is within 20% of aGM a mn/ (R a kg) (i.e. it scales with the gravitational poten- 
tial of the accreting WD before merger), with a = 0.270. Most of the equal and nearly 
equal mass mergers have convectively unstable cores. Artificially mixing these cores to 
make them isentropic decreases their maximum temperatures by 10-40% (see Fig. [I]). 

• The density at maximum temperature p(T max ) of the system is within 52% of 
central density of the donor, pd,c- 

• The central density of the remnant p c is within 67% of p ayC , the central density 
of the accretor. The central temperature T c for unequal mass mergers is significantly 
lower than T max . Since our conditions are far from the pyconuclear regime, temperature 
is more important than density, and so T max is more relevant to post merger evolution 
than T c , despite p c being higher than yo(T max ). 

• The enclosed mass inside the radius of highest temperature M enc (T max ) is 
within 35% of M a , the mass of the accretor, when Pd,c/p&,c < 0.6, i.e. for unequal mass 
mergers the hottest point is near the core-envelope interface. M enc (T max ) is zero when 
p&,c/ P&.c > 0.6, i.e. the hottest point is at the centre for nearly equal mass mergers. The 
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cut-off at Pd.c/Pa.c = 0.6 is sharp, suggesting that this ratio defines whether or not the 
merger is "equal mass" . 



4. Implications for the Sub-Mch Merger SN la Channel 

For Fig. [T] we have replicated |Van Kerkwij k et al.f s Fig. 1, and added our simulated 
merger remnants. Although several remnants are above the r cc = t v runaway line, none 
reach a carbon burning timescale shorter than the accretion timescale of several hours. 
Accretion of the disk onto the core is therefore of most immediate consequence. Accretion 
results in adiabatic compressional heating, and may increase density by up to an order 
of magnitude. Transposing the points on Fig. [T] one order of magnitude in density along 
a contour of constant entropy results in all remnants with accretor masses > 0.55 M 
reaching ignition (r cc = t v ). Those remnants with small mass ratios q < 0.5 will likely 
not explode, since as the runaway progresses, they would lose degeneracy around the 
time the fusion timescale becomes equal to the degeneracy timescale (r c i yn = r cc ). For all 
others, explosion seems, at least in principle, possible. 

A more accurate depiction of post merger evolution would require taking into consider- 
ation the the full viscous evolution, angular momentum transport, heating and compres- 
sion throughout the remnant (see Shen et al.|2011 for initial studies) . Likewise, the later, 



rapid nuclear runaway requires simulation to see if it can lead to detonation (Chang, 
White & Van Kerkwijk, in prep.). As noted earlier, merger remnants with their hottest 
point at the centre have Pd,c/Pa.c > 0.6 (corresponding to approximately M a — Md < 0.08 
M©). If this is required for an explosion, only a fraction of our remnants will explode. It 
is still possible, however, that some systems will experience off-centre detonations as a 
result of compressional heating. The heaviest of our equal-mass mergers may also explode 



as they merge, as described in Pakmor et al. (2010 2011 ) 



Nevertheless, Fig.[T]suggests that the sub-Mch merger channel for SNe la merits further 
investigation. 
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